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(N Abstract 



Repeated convolution of a probability measure on Z leads to the central limit the- 
orem and other limit theorems. This paper investigates what kinds of results remain 
without positivity. It reviews theorems due to Schoenberg, Greville, and Thomee which 
are motivated by applications to data smoothing (Schoenberg and Greville) and finite 
difference schemes (Thomee). Using Fourier transform arguments, we prove detailed 
decay bounds for convolution powers of finitely supported complex functions on Z. If 
M is an hermitian contraction, an estimate for the off-diagonal entries of the powers 
M£ of M k = I — (I — M) k is obtained. This generalizes the Carne-Varopoulos Markov 
chain estimate. 

> 

o 

^ 1 Introduction 

This is the first of a series of papers regarding the following question. What can be said 
about the convolution powers of a finitely supported function on a countable group? 

A lot is known if we are willing to assume that the function, say 4>, is non-negative and 
normalized by = 1. I n this case, the convolution power of is the probability 
. ^ distribution of the associated random walk after n-steps and estimating cp^ is a well studied 

^ problem. But what exactly is the role of "positivity" in classical results? 

This article focuses mostly on the simplest case, convolutions on the integer group Z. It 
reviews what is known (and why people worked on this question before) and provides some 
new results including an extension of a Markov chain estimate due to Carne and Varopoulos. 
In forthcoming papers, we will consider the case of Z d , which is significantly different from 
the 1-dimensional case, and the case of non-commutative groups where completely different 
techniques must be used and many open questions remain. 

Let 4> be a finitely supported probability measure on the integers Z. Define convolution 
powers of <\> by 
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The asymptotic "shape" of is described by the local central limit theorem (assuming 
irreducibility and aperiodicity): 



(n) (x) 



__ e -(x- an ) 2 /(2a 2 n) + Q 
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(1.1) 



In (1.1), o 1 = ^2,x 2 (j){x), a = ^x0(x), and the error term is uniform over x G Z. Proofs 
and many refinements are in EI] . 




Figure 1: 

For this and other basic theorems of probability, could there be similar results if (f>(x) 
is a more general function? For example, consider 

0(-2) = 0(2) = -i 0(1) = 0(-l) = ^ 0(0) = ^ (0(x) = otherwise). (1.2) 

The convolution powers 0( n ) are well-defined and results of Schoenberg |24j imply 

1 „ / x \ ( 1 



0W(a 



-ft 



(n/9)V4 ^VW9)V 



+ o 



1/4 



(1.3) 



with f/4 the real even function having Fourier transform e x . 

(-k -s-k 

More generally, H k is defined as having Fourier transform e - * for even k and 



for odd k. A graph of #4 appears in Figure 1 Also shown are the Gaussian density 



H<z(x) = 77j= e x 7/4 an d H$, H§. The function H3 is the famous Airy function (up to 



dilation) and it appears as the convolution limit of measures such as 

0(0) = 1 - 3a, 0(1) = 3a, 0(-l) = a, 0(2) = -a (0(x) = otherwise) (1.4) 

for < a < 1/4 fixed. Careful statements of general theorems are in Section 2} We note 
that there are fundamental differences between the functions and the function f^fc+i- 
An exact formula exists only for H2 but H2k is even with ff2fc(0) > and satisfies 

Vx, \H 2k [x)\ <C k eM-c k \x\ 2k/{2k - l) ). (1.5) 
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In particular, H 2 k is in L 1 (M). The functions H 2 k+i satisfy i?2fe+i(0) > but are not 
even, decay differently along the negative and positive semi-axes and are not absolutely 
integrable. For any k > 2, changes sign infinitely many times. 

The following theorem illustrates the main result of this paper. It gives an exponential 
upper bound of the type ( 1.5 ) for the convolution powers of certain complex valued functions 
on Z. Such bounds are implicit in the literature on finite difference methods. See [28\ I29j. 
The second part of the theorem gives lower bounds on the real and imaginary parts in an 
appropriate neighborood of 0. 

Theorem 1.1. Let<j)(x) be a finitely supported complex-valued function on7L with^2 x (j)(x) = 



1. Assume that (f)(6) = ^ a 



[x e 



1X0 



satisfies 



<1 for 6 € (— 7r,7r], 6»/0. 



(1.6) 



Assume further that there exist an even integer u and a complex number 7 with Re (7) > 
such that 

0(0) = e -^d+o(D) as 9^0. (1.7) 
Then, there are constants c,C£ (0, 00) such that, for all x £ Z, n £ N*, 



c/> (n) (a 



C 

< —^-r- exp 



n 



X/u 



x 



n 



l/v 



Further, there are constants c±,C2 > such that, for alli£2,n£N with \x\ < c\n l / u , we 
have 

Re(^ (n) (x)) > c 2 n- 1 / v 



and 



• «/Im( 7 )^0, 

\lax((j)(x))\ > c 2 n- l ' v 
with \m.(cj)( n \x)) having the same sign as 1111(7) ; 

• if Im(7) = 0, 

\lm(<t>{x)\) = o(n- x ' u ). 



The first part of Theorem |1.1| is proved in Section 3 The second part follows from a local 
limit theorem derived in ISection 2.31 This local theorem gives further information about 



(j)( n ) i n the region |x| < c\n x l v . Each of these results is obtained using Fourier transform 
techniques The conclusions of Theorem |1.1| can be described roughly as follows. Under the 
given hypothesis, |0( n )(2;)| attains a maximum of order l/n 1 /^ and is bounded below by 
l/n l / u in a n^-neighborhood of 0. Further, as x moves away from 0, |<p n )(a;)| presents a 
relatively fast decay controlled by 



exp 



\x\ 
n x l v 



In modern random walk theory, when <p is a centered probability and v can only take 
the value v = 2, the behaviors described in Theorem |1.1| are often discussed under the 
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names of "diagonal behavior" (more precisely, "near diagonal behavior") and "off-diagonal 
behavior" . In particular, in the theory of random walks on graphs, long range off-diagonal 
upper bounds of the type e~ cd '( x ' y ' ' n are known as Carne-Varopoulos bounds (here d(x, y) 
denotes the natural graph distance). In Section Section 4 we explain how Theorem 



leads to generalizations of Carne-Varopoulos bounds. 

In general, it is not true that the £i-norm of the convolution powers <f>^ of a complex 
valued function (f> with ^0 = 1 stay bounded uniformly in n. However, as an immediate 
corollary of Theorem [LTj we recover the following known result. 



Corollary 1.2. Under the hypotheses of Theorem 1.1 there exists a constant C such that 



Vn, ^|0(")(x)| < C < oo. 

This important property is studied in detail in [IJ [281 USB EE] because of its close connec- 
tions to the "stability" of certain approximation schemes for partial differential equations. 
See Theorem 12.11 below and the short discussion in lSection 5.21 

The hypotheses made in Theorem |1.1| are somewhat mysterious. Indeed, it is very 
unclear how to replace these hypotheses in the context of non-commutative countable groups 
(for which no viable Fourier transform exists, in general). Even on Z, these hypotheses are 
not entirely natural since the Fourier transform of a function (j> with Y2 X 4>(x) = 1 may 
well attain its maximum at multiple points and not at 0. Sections 2 and 3 discuss this in 
detail. 

We end this introduction with a brief description of the content of the paper. Section 2 
is devoted to local limit theorems for convolution of complex valued functions. In particular, 
Section 2.3| provides an extension of earlier results of Schoenberg and Greville and |Section 2.4 



describes illustrative examples. Section 3 is the main section of this paper and treats upper 



bounds of the type (1.8). See Theorems 3.1 and 3.4. In Section 4, the results of Section 3 
and the discrete transmutation formula of Carne are used to obtain long range estimates for 
the powers of the operator = I — (I — M) k when M is the infinite matrix of a reversible 
Markov chain (more generally, an hermitian contraction). Section 5 gives pointers to various 
earlier works where convolutions of signed measures play an explicit or implicit part. 



2 The results of Shoenberg, Greville and Thome 

In this short section, we briefly review results by Schoenberg, Greville and Thomee. Schoen- 
berg and Greville where motivated by the earlier work of de Forest and statistical data 
smoothing procedures. Thomee's work is motivated by numerical approximation schemes 
for differential equations. Brief explanations regarding the motivations of these authors and 
other applications are collected in Section 5. 

2.1 Stability 

The largest body of work concerning the problem discussed in this paper can be found in 
the literature concerning numerical approximation schemes for differential equations with 
work by John |18| . Aronson [1], and Widlund |31| . Thome's articles |28t 129] give excellent 
pointers to the literature. However, these papers do not isolate or emphasize the basic 
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convolution aspect of the results. This makes extracting the relevant results somewhat 
difficult. Nevertheless, the following essential result is explicitly stated in [29]. Given a 
complex valued absolutely summable function on Z, consider the property 



Vn, ^|0^(x)| < C < oo. (2.1) 



Theorem 2.1 (See |29t Theorem 7.54]). Condition (2.1) is satisfied if and only one or the 
other of the following two conditions is satisfied: 

1. 0(£) = (e^ for some y G Z and ( G C twtfc |C| = 1 

,2. |0(C)| < 1 except for at most a finite number of points £ g> q G {1, . . . , Q} in |£| < 7T 
where |a(f g )| = 1, and there are constants a q ,fi q ,v q , q = 1,...,Q, with a q real, 
Re(Pq) > and u q an even natural integer, such that 

kt q + = ktq) exp (ia q i - p q C« (1 + o(l))) as £ -> 0. 



Thomee also describes what happens when the principal nonlinear term of the expansion 
at a point is imaginary (of odd or even degree). In particular, if |0| < 1, 0(0) = 1 and is 
the unique point where |0| is maximum with 

0(0 = exp (ia£ + ij^ - PC (1 + o(l))) as 
with a, 7 real, 7 ^ 0, Re(/3) > 0, < fj, < v and v an even integer, [291 (7-9)] gives 

l^ (n) (f )l - nC 1 -"/")/ 2 . (2.2) 

y 

Here, ~ means that the ratio of the two sides stays between two positive constants as n 
tends to infinity. 

Condition (2.1) is called a stability condition. In the case of smoothing, it shows that 
the iterated smoother is continuous from £°° to £°° ; small changes in the input sequence lead 
to close smoothed sequences, uniformly in the iteration number n. In the case of divided 
difference schemes, it shows that the scheme applied to bounded data gives bounded output, 
uniformly in time. 



2.2 De Forest local limit theorems 

Following de Forest, Schoenberg [24J and Greville |16j explore local limit theorems for real 
valued that are allowed to change sign and are normalized by ^ = 1. Schoenberg treats 
the case of symmetric (i.e., 0(x) = 0(— x)) under the assumption that |0(6>)| < 1 for 
9 G (0,2vr) and 0(0) = 1 - X9 k + 0(|6>| 2fc+1 ) (with, by necessity, A > 0). Greville observes 
that symmetry is not essential and also treats the case when0(6») = l + ai0 2k+1 + 0{\6\ 2k+2 ). 

For each positive integer k, let denote the function of one real variable defined by 
Hk(0 = e - ^ for even k and Hk(0 = e 1 ^ for odd k. 

Theorem 2.2 (Greville |16) and Schoenberg [24], in the spirit of de Forest). Let (ft be a real 
summable function on Z with Fourier transform satisfying \4>(0)\ < 1 for all 9 G (0, 2ir) 
and 

0(0) = 1 + a(iB) m + O(|0| m+1 ) with 0/a£l and m an integer greater than 1. (2.3) 
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• If m is even then —ai m = A must be positive and we have 

= (An)" 1 /" 1 ^ (xiXn)- 1 '™} +o( n - 1/m ), 

uniformly over x € Z. 

• If m is odd and e denotes the sign of the real ai m ~ l , e = sign(az m_1 ) ; then we have 

0M( X ) = (\ a \n)~ 1/m H m (ex(|a|n)" 1/m ) + o(n" 1/,n ). 
Remark. If we consider the expansion 4>{9) = ^ ■ Aj(i9y then A,- = i ^ x (j){x)x 3 . Hypoth- 



esis (2.3 ) is simply the assumption that the first k— 1 moments of <f> vanish and a = ^ 0. 



Note that a is real here because <f> is real. 

Remark. The hypotheses made in Theorem |2.2| are somewhat ad hoc. Schoenberg and 
Greville are interested in data smoothing procedures. This explains the basic hypotheses 
that cp is real and satisfies ^ <f> = 1 but these properties are not essential, at least in the 
case when m is even. The assumption that max|^| = 1 is actually the essential condition. 
Coupled with the condition ^ <p = 1, it implies that \<j)\ attains its maximum at 9 = and 
that, in fact, 0(0) = 1. Schoenberg and Greville then make the additional assumptions that 



< 1 if 9 ^ 0. The condition (2.3) simply captures the vanishing order of <f> — 1 at 0. 
It excludes the possibility of a drift. In the next section, we extend the even m result to 
allow complex (J), a drift, and multiple points where \4>\ attains it maximum. 



Remark. The proof of Theorem 2.2 is significantly more difficult in the case where m is odd 



than in the case when m is even. Indeed, for odd m = 2k + 1 > 1, 

1 f°° 

H 2 k+i(x) = — cos(n 2fc+1 — xu) du. 



o 



This integral is not absolutely convergent and i?2jfc+l is neither integrable nor square inte- 
grable. Nevertheless, it is an entire function. 

2.3 An extended local limit theorem in the even case 



Theorem 2.2 is taken from [16]. Its form is somewhat constrained, perhaps due to the 
applications that Greville had in mind. Comparison with Theorem |2.1| leads to the fol- 
lowing extension where H2k,b denotes the complex valued Schwartz function whose Fourier 
transform is H 2 k,b(0 = e~^ 1+ib ^ 2k , b € R. 

Theorem 2.3. Let <p be a complex absolutely summable function on Z with Fourier trans- 
form 4> satisfying \4>(0)\ < 1 for all 9 £ [— vr,7r]. Assume that \4>(£)\ < 1 except for at most a 
finite number of points £ 9 , q £ {1, . . . , Q}, in |£| < ir where \4>(£,q)\ = 1- Assume further that 
there is an integer Q\ < Q, an even integer m and constants a q , j3 q = b q +ib', q = 1, . . . , Q\, 
with a q , b q , b' G M, b q = Re(f3 q ) > such that 

V<?e{l,...,Qi}, 0(6 z + £) = <toexp(m g £-A ? r (l + o(l))) as £ ^ 0; (2.4) 
and 

Vg€{Qi + l,...,Q}, m q + 0\ <exp(-c g |en as £ -+ 0, 
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for some p G (0, m). Then we have 
Qi 

(x) = Y,{\n)- l l m e-^^ q ) n H mK/bq {{x - a q n){b q ny l l m ) + o(n^ m ), 



1 

where the error term is uniform in x G Z. 



Remark. To make sense of this result, it is necessary to know something about the functions 
H 2 k,/3- By definition, 

1 f°° 

Re(H 2kfi (x)) = - (cosxO(cos(^ 2fc ))e-« d£ 
K Jo 

and 

lm(H 2k , b {x)) = - / (cos^(sin(&£ 2fc ))e-^. 

From these formula, one easily sees that Re(H 2 k,b(ty) is always positive and that Im(i7 2 fc,b(0)) 
is zero, positive or negative if and only if b is zero, positive or negative. Further, it holds 
that 

2k 

\H 2 k,b(x)\ < Cfc,&exp(-c fcjfe |x| 2fc -i)- 

In the case k = 1, we have 

1 M 2 
H 2b (x) = , e 4 ( i+ib). 

This complex valued function of the real variable x is the heat kernel, i.e., the kernel of e zA , 
computed at the complex time z = 1 + ib. Here A denotes the unique self-adjoint extension 
of (d/dx) 2 originally defined on smooth compactly supported functions. 

Proof. Set / = (— 7r,7r]. For each q, let I q = [i q — e, £ g + e], e > 0, be a small interval centered 
around £ 9 and let J = I\ufl q . We have 



27r<A (n) (x) = e- ix9 [0{6)] 

J — TV 

= I e- ix9 [(j)(9)} n d9 + Y^ ( e- ixe [j>{9)] n d9 = J2® 

J J <7=1 



*d0 



where 



$ = /" e -" e [^)] n ^ and / e _ia *[0(0)] n d0, g= 1,...,Q. 

J J Jl a 



On J, there exists p = p e G (0,1) such that |0| < p and thus, |3>o| = 0(p n ). For q G 
{Qi + 1, . . . , Q}, we have 



\$ q \ = | y" e-^'+^fo + fl)]"^ 
< J 6 e~ c ^d9 

/oo 
e - c «H M du = o(n- 1 / m ). 
-oo 
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The main contribution comes from the integrals & q , q £ {1, • • • ,Qi}- We set 

^,,n(«) = V ^") -1 '™^, + u(b q n)- l / m )] n , 



Uq,n 



x — a„n 



(b q n) l l m ' 



and write 



(b q n) 



(fegn) 1 /"' 

■(6,11)!/" 



e- iy ^ u ^ n (u)du. 



Next, since 



we have 



^H m y/ bq (y q)n ) 



-y q ,nu e -{i+ib' q /b q )u m ^ u 



% - 27r(b q ny 1 / m e-^fa q ) n H m , Vbq (y q , n )) 



f e- iy «'" u ij q , n ( U )du 



-Wq,nU e -(l+ib' q /b q )u m ^ 



e- iy "- nU [^ n (u) - e-( 1+ib '«/ b ^ u "']du 

e -iy«,nV e -(l+iW b <l) um du=ll -1 2 . 



\u\<e{b q ny/ r - 



, \u\>e(b q n) 1 / r ' 

The integral X 2 can be estimated brutally by 



\X 2 \ <2 f e- um du = 0(e- emb « n ). 

Ju>t(b q n) 1 / m 



To estimate Xi, note that (2.4) shows that for any rj > there exists e > such that 
for all |u| < e(6 f/ n) 1 / m . It follows that 

/•oo 

|Xa| < / e-^/^u. 

JO 

Putting all these estimates together, we find that 

<P^(x) - (6 ? n)- 1 /- e -^0(^)^ mj6 , /b9 ( %n )) = (l/ n Vm). 



□ 
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2.4 Examples 

One of the simplest example that can be used to illustrate the results discussed in this paper 
appears in the following proposition. 

Proposition 2.4. Assume that <f> is real symmetric, \<j>\ < 1 and there exists a > such 
that 

0(9) = l-a6 2k {l + o(l)) at 6 = 0. (2.5) 

Then <p(x) ^ for some x > k. If we assume further that <f> is supported on {—k, . . . , k} 
then a G (0, 2~ 2fc+1 ] and 

4> = S - X(S - (3) k = S - (S - (3 xl/k ) k (2.6) 

= B-i)^(-)^ 

where A = a2 k G (0,2- fe+1 ] ; (3 = §($_i + 6{) and (3 s = (l- s)S + sf3, sG (0,1). 
Remark. The function <f> defined at (2.6) satisfies 

fad) = 1 - A(l -cos9) k . 



It follows that max{|</>|} = 1 as well as (2.5) for any A = a2 k G (0, 2 _fc+1 ].The maximum 
1 = max{|0|} is attained solely at if and only if A G (0, 2~ fc+1 ). Note that any parameter 
A in the range (0, 1/2] is admissible for all values of k. 

Remark. If convolution by So — (3 S (for some fixed s G (0, 1), say s = 1/2) is interpreted as 
the discrete analog of the (positive) Laplacian — d% then convolution by 5q — <j> = (So — Ps) k 
is analogous to the higher even powers of the Laplacian, that is, (— l) k d 2k . 

Proof. By assumption, if p is the largest integer such that <j)(p) ^ 0, <f> is a polynomial Q in 
cos# of degree p. If (2.5) holds then the polynomial 1 — Q vanishes of order k at 1. Hence, 
we must have p > k. If we assume that <f> is supported on {—k, . . . , k}, then we must must 
have 1 — Q(u) = A(l — u) k and a = 2 k X with a as in (2.5). That is, 



0(0) = 1 - A(l - cos#) fc = 1 - (l - (l - A 1/fe ) + X 1/k 



cos ( 



The condition that \(f)\ < 1 on (0, 2n) translates into 1 - X2 k > -1, that is A < 2~ k+l . Fur- 
ther, since cos 9 = (3(9) where j3(l) = f3(— 1) = 1/2 and (3 = otherwise (i.e., /3 =Bernoulli(l/2)), 
we see that we must have 



i=i v 7 



where f3 s = (1 — s)Sq + sj3. □ 



Example 2.1. In Proposition 2.4, consider the case whe n k = 1 and A = 1 so that <f> = j3 



We have (f>(9) = cos 9. In this classical case, Theorem 2.3 yields 

^)( x ) = (1 + (-l)-e— )(n/2)- l l 2 H 2 (x/(n/2) l l 2 ) + o^ 1 / 2 ). 
This captures the periodicity of the Bernoulli walk. 
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I I (2) 

Example 2.2. In Proposition 2.4 consider the case when k = 2, t hat is = 2(3 S — fis 



with s G (0, 1/V2). We have 0(0) = 1 - s 2 (l - cosfl) 2 . Theorem \2.2\ yields 

(n) (x) ~ (4/s 2 n) 1 / 4 J ff 4 (x(4/(s 2 n) 1/4 ). 



The same result holds true in the limit case where s = l/v2 and = 2/3-^^ — Z^//^- 

However, in this case, |0| < 1 and (f>(0) = 1 if and only if 9 = or 9 = tt. At 0, 0(0) = 
1 - ||0| 4 + 0(\9\ 6 ). At tt, 0(0) = -1 + (tt - 0) 2 + 0(|vr - 0| 4 ). To obtain the desired 
asymptotic, apply Theorem |2.3| Compare to the previous example. 



Example 2.3. Let be defined by 0(0) = 5/8, 0(±2) = -1/4, 0(±4) = -1/16 and 
0(z) = otherwise. We have 

4,(9) = I - Jcos2#- ^cos4# 
8 2 8 

5 1 

= - - cos 2 9 - - cos 2 26* 
4 4 

= 1 - cos 4 9. 

Hence |0| < 1 and 0(0) = 1 if and only if 6» = ±tt/2. Further, at 0± = ±vr/2, 

0(0) = 1- \9-9±\ 4 + 0(\9-9±f). 



Theorem 2.3 applies and gives 

0(™)( x ) = ( e -**f +e ix ^)n- 1 ^H A (xn- 1 ^) + o^ 1 / 4 ) 
= 2 cos(vrx/2)n- 1 / 4 J ff 4 (a;n- 1 / 4 ) + o^ 1 / 4 ). 

Example 2.4. In this example, we consider the convolution powers of an arbitrary real 
valued function supported on {—1, 0, +1} (except for some trivial cases). For a®, a + , a_ € E, 
let be given by 

0(0) = aoj 0(=tl) = a ± and = otherwise. 

We assume ao > and that either a+ 7^ or a_ 7^ to avoid trivialities. We do not assume 
any normalization. In particular, 

^ = a + a + + a_ 

is an arbitrary real number. The following proposition shows that there are essentially 3 dif- 
ferent "generic" cases (each occurring on an open subset of the parameter space) . In only one 
of these 3 cases is the normalization ^0 = 1 the correct normalization giving max{|0|} = 1. 
In the other two cases, different normalizations are needed to insure max{|0|} = 1. Because 
we do not incorporate any normalization, the asymptotic described below for 0( n ) contain 
an exponential term A n . In each case, the constant A satisfies A = | max{|0|} and is given 
explicitly. 

Proposition 2.5. Referring to the function defined above, the following asymptotics hold 
true: 
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Assume that a + a_ > or that a + a_ < and 4|a + a_| < ao\a + + a_|. S'ei 

a+ — a_ |a + | + |a_| a 2 
A = a + |a+| + |a_|, a = -r and 7 = — — . 

T/ien 

0(n)( x ) = f ^-) 1 A»( 7 r 1 )- 1 /2 e -l-H 2 /7« + o(iV 1 / 2 ) 
\|a+ + a_| / 

where the error term is uniform in x. 

Assume that a+a_ < and 4|a+a_| > ao|a+ + a_ |. Se£ 

A= |a+ -a_|(l + ag/4|a+a_|) 1/2 , a + "~ 



a+ — a_ 



£e£ #0 £ (0) 71 ") ^ e defined by cos#o = — «o( a + + a_)/4a + a_ and set 

4ao|a + a_| sin#o 16|a + a_| 2 — On(a + + a-) 2 

(a+ - a_)(4|a + a_| + a 2 ,)' 2(ag + 4|a+a_|)(a+ - a_) 2 ' 

Zei wo 6e i/te argument of 4>{0q). Then 

^ n \x) = (n/ 1 )- 1 / 2 A n e- ixdo+inuJo H 2tbh ((x-an)/(-rn) 1 / 2 ) 

+(n/ 1 )- 1 / 2 A n e ixd °- inuJ °H 2 ,_ bh ((x - an)/( 7 n) 1 / 2 ) + o{A n n - 1 / 2 ) 

where the error term is uniform in x. 

Remark. The two principal terms on the right-hand side of the last equation are complex 
conjugate so that their sum is real (as it should be since <fi is real valued). 

Remark. This example can also be used to illustrate the stability theorem, Theorem |2.1| 
Indeed, all cases with qq > and either a+ or a_ non-zero are considered in Proposition 



2.5 except for the very special case when 4|a + a_| = ao|a + + a_|. Let A = A(ao, a+, a_ 



be as defined in Proposition |2.5| As a corollary of the proof given below, it follows that 



the normalized function 0o = <fi/A satisfies the stability condition (2.1) in all cases but the 
special case 4|a + a_| = ao| a + + a -| f° r which it actually fails. 

Proof. Obviously 

4>{6) = ao + (a+ + a J) cosO + i(a+ — a-) sm# 

and 

\4>(9)\ 2 = al + a 2 + + a 2 + 2a + a_(2cos 2 9 - 1) + 2a (a+ + a_)cos6». 

Further 

0(0) = ao + a + + a_, 0(7r) = ao — (a+ + a_). 

If a+a_ > then \<j)\ has a maximum which is attained only at if a + + a_ > and only 
at 7r if a + + a_ < 0. 

Assume first that a + a_ > and set 

a + — a_ |a + | + |a_| (a + — a_) 2 
A = a + |a+| + |a_|, a = and 7 = — — % > 0. 
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Considering separately the two cases a + + a_ > and a + + a_ < 0, we obtain 

,2 



e-^UiO)' 1 ^) = {l-ia6-^-9 2 + o{\9\ 2 )){l + ia9-{ 1 + a 2 /2)e 2 + o{\9\ 2 )) 

= l- 7 ^ (1 + o(1)) = e -7^(l+o(l)) > 

Hence, if a + + a_ > we have 

<^ n )(x) = A^n)" 1 / 2 ^^ - an)/( 7 n) 1 / 2 ) + o^n" 1 / 2 ). 

If instead a + + a_ < then we have 

= {-lfA n { 1 n)- 1 ' 2 H 2 {{x - an)/^) 1 / 2 ) +o(A n n" 1 / 2 ). 

Next we consider what happens when a + a_ < 0. Computing the derivative of f(9) = 
1 4>(6)\ 2 gives 

f{9) = -2(4a + a_ cos# + a (a + + a_)) sin0. 

If 4|a + a_| < ao|a+ + a-|, then |0| 2 attains its maxima at if a + + a_ > and at 7r if 
a + + a_ < 0. In each case, the asymptotic is the same as described above. 
If 4|a+a_| > ao\a+ + a_|, then we set #o to be the point in (0, ir) such that 

o (a+ + a_) 

cos do = - 



4a+a_ 

and note that \4>\ 2 has twin maxima at 9 = ±#o where 

A 2 = |^(±g )| 2 = a l + a 2 + a 2 + 2[Q+a _ | + Q 2 (g+ + 

4| cz-)-£i — | 

= (a+-a_) 2 fl + °° 



4|a + a_ 



Further 



Re(<^(0o)0(0)) = (ao + (a+ + a_) cos0o)(ao + (a+ + «-) cos0) + (a + — a_) sin0osin6* 



a 2 

(a+ — a_) 2 I — — ^ — . + cos 9q cos 9 + sin^o sin( 
\4\a + a- 1 

(a+ - a_) 2 f - °° . + cos(0 o - 0) 

^'K'-^b' 9 - 90)2(1 
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and 

lm((/>(0o)4>(9)) = (a+ — a_)((oo + (a+ + a_) cos#o) sin# — (ao + (a+ + a-) cos 6*) sin#o) 

= ao(a+ — a_)(sin0 — sin#o) + (a+ — a?-) sin(# — #o) 

= (a_|_ — a_)(oo cos #o + «+ + a -){® ~ $o) 

sin0 o (#-#o) 2 (l + o(l)) 



o (a+ - a_) /( in n i2 



2 



■sinfl (#-#o) 2 (l + o(l)) 



«o(Q+ ~ a -) „ ; ,. ,„ ,2, 

2 



a.f — o_ 
|2 4ao|a+a„| sin^o 
(o+ - a_)(4|a + a_| + a 2 ,) 



IM»)l a ,_ T'^T-'T.. a , (g - ^o) 2 (i + o(i)) 



Set 



and 



a + + a_ 4ao|a_|_a__| sin^o 

a = > 0=7 7771 1 97 

a + — a_ (a + — a_J(4 a + a_ + %) 



2|a + a_| 1 / a + + a_ x 



aQ + 4|a+a_| 2 ^0+ — a_ 
16|a + a_| 2 — ao(a+ + a_) 2 
2(a 2 + 4|a+a_|)(a+-a_) 2 ' 

Note that 7 is (strictly) positive if 4|a+o_| > ao| a + + a —\- With this notation, assuming 
that 4|a + a_| > ao|°+ + a -|> we have 

o + 0) = ^ ) e ^-(7+^ 2 (i+o(i)). 



Similarly, 

£> »a0-(7-i/8)0 2 (l+o(l)) 



+ = -00 e" 



2.3 



gives 



Let wo be the argument of 4>(9q). Then Theorem 

</,(») (x) = (n/7)- 1 / 2 A n e-^ 0+inwo ^ 2iV7 ((x-an)/( 7 n) 1 / 2 ) 

+ (n/ 7 )- 1 / 2 ^e^- inwo F 2i _ 6/7 ((x - an)/( 7 n) 1 / 2 ) + o^- 1 / 2 ). 

Finally, if 4a + a_ = — ao(a + + a_) (resp. 4a + a_ = ao(a + + a_)), |0| is maximum 
at (resp. at tt). The two case are similar and we treat only the case when 4a+a_ = 
— ao(a + + a ). In this case we have 

^(O)- 1 ^) = e ia0-i|(a-a 3 )0 3 -|(a 2 -a*)^(l+ O (l)). 

□ 

By J2T2I, in this case we have £ \<fi n \ x )\ ~ A^n 1 " 3 / 4 . 
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3 Bounds on convolution powers of normalized complex func- 
tions on Z 



The goal of this section is to give good upper bounds for the convolution powers (ft^ of 
a given complex valued function that is finitely supported on Z. Let (ft be the Fourier 
transform of (ft so that 

<hW(x) = ^- [ e- ixe \<b(0)]",W. 



2tt 

Obviously the function (ft can be normalized in some appropriate way and it is very reason- 
able to chose the normalization 

max{|0|) = 1. 
eel 

Note that this is the same as saying that the operator norm of the convolution operator by 
(ft acting on £ 2 (Z) is 1. 

As we assume that (ft is finitely supported, either = 1 on /, that is, (f)(9) = ce ij with 
|c| = 1, j an integer, or max{|0|} is attained at at most finitely many points in / (the zeros 
of the trigonometric polynomial \(j)\ 2 — 1). 

In the first subsection of this section, we concentrate on the case when max \(ft\ is attained 
at only one point 9q £ I and (f)(0) — 4>(9q) vanishes up to some even order. 

The case when the maximum of (ft is attained at more than one point will be considered 
in the second subsection below but our results are much less precise in that case. 



3.1 Generalized Gaussian bounds, I 

In this section, we assume that max|0| is attained at only one point 9q £ I- By replacing 
(ft by x ■-> (j) 9o (x) = ^(6»o) _1 e i:Eeo 0(x), it is enough to consider the case when 9q = and 
(f)(0) = 1. So assume from now on that (f> is finitely supported with 

{x : 4>(x) 7^ 0} C [—K,K], Wi = 2M«» 

and 

0(0) = 1, |0(0)| <1 for 9 el* = J\{0}. (3.1) 

Since (ft is finitely supported, (ft is actually an entire function of the complex variable z = 
u + iv E C, namely, 

^(z) = ^Hx)e izx . 

As noticed in [28, 29] (and many other places), it is most efficient to expand log (ft and write 

oo 

(ft(9) = exp(^^Cj# J ) near 

i=i 

to study the convolution powers (ft( n \ Indeed, the conditions 0(0) = 1 and \(ft\ < 1 in I* 
implies that only two cases may arise. Namely, either there exist a real a, two integers 
< /j, < v, v even, a real polynomial q with q(0) ^ 0, and a complex number 7 with 
Re (7) > such that 

£fQ\ _ e ia0+;0^(<?)-70"(l+ O (l)) / 3 _ 2 ) 
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or there exists a real a, an even integer v and a complex number 7 with Re(7) > such 
that 



o ia0-70" (l+o(l)) 



(3.3) 



When (3.2) holds, [2H] shows that ||<^ n )||i tends to infinity with n and in fact 



See |29l Sect. 7] and the references therein. Hence we focus on the case when (3.3) holds. 
The simplest form of our main result is stated in the following theorem. 

Theorem 3.1. Let cf) : Z -+ C be a finitely supported complex valued function such that 
(3.3) holds true and \4>{0)\ < 1 on I*. Then there are constants C,c£ (0, 00) such that for 



any x 6 Z and n £ N* ; we have 



n 



C 

"TA7 



exp 



\x — an\ 



n 



l/v 



Proof. Write 



where 



2tt 



-ixO 2 ( a\n 



<t){e) n de 



2tt 



e- l{x - an)e P{6) n d6 



P(9) 



-iad 



Condition (3.3) together with \4>\ < 1 on I* implies that 

P(0) = exp(-7^(l + o(l))), \P\ < 1 on I* 



(3.4) 



with 7, v as in (3.3), that is Re(7) > and v is an even integer. In addition P is obviously an 



entire function of the complex variable z = u + iv and (3.4) implies that there are constants 
7i;72 £ (0,oo) such that 



|-P(-z)| < exp(— 71?/' + 72f") on {z = u + iv : \u\ < 37r/2} 



and 



(3.5) 
(3.6) 

To see why (3.6) holds, at the origin, use the assumed expansion P(u) = e * Q «-7« I '( 1 + ( 1 )) - 
Away from the origin, note that 4>(z) is a trigonometric polynomial and thus is periodic in u 



\P(z)\ — exp(72f^) on C. 



with growth at most exponential in \v\ for large v. The estimate (3.5) improves significantly 



upon (3.6) only in the sector 



{z = u + iv : \v\ < —\u\; \u\ < 37r/2} 
7i 



where the the two parts of the hypothesis (3.4) together give the desired result. 
Next, write 



^ n \x) 



1 



■ (x-an) , 



n 1 /" P n {6) d0 
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with 

It follows that 



n 



\x — cm 



V"j = (-n^vr, n^vr] and P n {9) = P{9/n l l u ) n . 



n 



l/v 



\^\x)\ 



27m 1 / 1 



e n 1 /- u d q e P n (9) dO 



(3.7) 



i (x-an) g 

This step requires some explanation because when an is not an integer, neither e i 1 /" 
nor P n (9) are periodic function of period 2ixn 1 l v so that, a priori, the repeated integration 



by parts used to obtain (3.7) should produce boundary terms. However, introducing the 



operator 5 defined on smooth functions by 

5f(9) = e- l ^ e d e [e^ e f(9)], 

write 



. (x-an) , 



n 1 /"! 



[die °}P n (9) d9 



e « 



an a ■ an n ■ x n _. , 

d q Je ^/" e ni/»°](<p(9/n l / u )) n d9 



5 r [e i ^ e ](^(9/n 1/u )) n d9. 



Next, observe that 6 preserves the periodicity of the function /. In particular, if / is 
periodic of period 2im 1 l 1 ' then 5f has the same property. Hence the formal adjoint of 5 on 
27rn 1//iy -periodic functions is given by 



5*f(9) = -e^»d e [e V/*7(0)]. 



It follows that 



5 r [e V/- e ]($(9/n 1 / v )) n dO = {-If e ^ V^fl/n 1 ^))™ d9 

= (-l) r / e- 1 ^^ 9 d q e P n {9) d9. 



This justifies (3.7). We note that a formula similar to (3.7) appears on page 124, equation 



(2.2), of the classic paper by Ney and Spitzer [20J without warning or detailed explanations. 
Observe that (3.5)-(3.6) translate immediately into 

\P n (z)\ < exp(-7iu 1/ + j 2 v v ) on {z = u + iv : \u\ < 3n 1/v vr/2} (3.8) 

and 

\P n (z)\ < exp( 72 ^) on C. (3.9) 
The crucial estimate is given by the following proposition. 

Proposition 3.2. There are constants C\,a\ S (0, oo) such that for any q,n = 0, 1, 2 ... , 

and 9 € n x l u I, 

\d q e P n {0)\ < Cl+iqlq-^expi-a^)- 
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Assuming that this proposition has been proved, we obtain 
\x — an I 



n l/v 

This is of the form 

\0 {n) (x)\ < C 2 nr 1/v M- g q^- 1 ^ q , M 



\x — an 



Elementary calculus shows that 

inf {M' q q^~ l/v)q } < Cexp(-cM^r), c 



g =o;i; 2 ,... L * j-— ^ „ e{ v -iy 



This gives the upper bound stated in Theorem |3.1[ 

It remains to prove Proposition |3.2| By Cauchy's formula, 



'151 

Consider two cases. 

If q < 9 U , pick r = eq l / u with e > small enough (depending on 71,72 in (3.8)) so that 



|P n (z)| < exp(— jsO u ) on \z — 9\ = r. This easily gives the inequality of Proposition 3.2 
when q < 6 U . 

If, instead, q > U then pick r = q l / u and observe that, on \z — 6\ = r, 
\Pn{z)\ < exp(27 2 r y ) < exp(-ai0" + (2 72 + ai )q). 
This yields the desired estimate when q > 6 U . □ 



A useful complement to Theorem 3.1 involves "regularity" estimates for cf)( n \x). Namely, 
for any integer y and function / : Z — > C, set d y f(x) = f(x + y) — f(x). 



Theorem 3.3. Let <f> '■ Z — > C be a finitely supported function such that (3.3) holds true 
and \4>(9)\ < 1 on I*. Then there are constants A,C,c G (0, 00) such that for any x G Z, 
n G N* and any y\, . . . , y m G Z with \yj\ < An 1 ' 11 , j = 1, . . . , m, we have 

\d ■■■d d>W(x)\ < <Jm MM exp ( ( \x-c*n\ 
\°yi °ym<P { x )\^ n (i+m)/u exp I c y n \jv 

In particular, there exists C\ such that for all n > 1 and y±, . . . , y m G Z with \yj\ < An 1 ^ , 
j = 1, . . . , m, we have 
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Proof. Observe that 



1 



2irn 1 / u J n i/vi 



■ (x-an) , 



e l ^r u Q n {9)de 



where 



3=1 



Clearly, for e > small enough, if z = u + iv with \u\ < 3n 1 ^7r/2 and \v\ < e\u\, (3.5) gives 



K3 =1 



Further for all z G C, (3.6) gives 



u' =1 



Using these estimates in place of (3.5 )-( 3.6 ) and the line of reasoning of the proof of Theorem 
3.1 gives the desired estimate. □ 

We end by stating explicitly the most general result obtained in this section. 

Theorem 3.4. Let <f> : Z — > C be a finitely supported function such that there exists 6q G 
/ = (— 7T, 7r] such that 

\$(6 )\ = 1 and \$(6)\ < 1 on I \{6 }. 

Assume further that there exist a £ R, 7 £ C with Re{^) > 0, and an even integer v such 
that 



<i>(9a + 0) = <f>{e )e 



iaG-~i6 v {l+o(l)) 



at 9 = 0. 



(3.10) 



Then there are constants A,C,c G (0, 00) such that for any x G Z, n G N* , m = 0,1, 
and yi , ... , y m G Z with \yj\ < An l / U , j = 1, . . . , m, we have 



)- n e ixdo ^ n) {x)\\ < 



c m n?\y3-\ 



exp — c 



n (l+m)/u 

Further, there exists e, rj > such that, on \x — an\ < en l / v , 

Re (cK# rV*V n) (x)) > r)n- l l v . 



\x — an\ 



n 



1/v 



Proof. This follows from Theorems 3.l||3.3 applied to (pooi 3 
assertion follows from Theorem [2JJJ 



' )- L e lx0 °4>(x). The last 

□ 
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3.2 Sub-Gaussian bounds 



Examples such as 2.4 show it is very natural to allow the Fourier transform 4> to attain its 
maximum of 1 at more than one point on (— vr, ir]. This section explores briefly what global 
estimates can be obtained in such cases. In the following theorem, <j) is allowed to attain its 
maximum, 1, at finitely many points. At each of this points, we assume that the dominant 
non-linear term in the expansion of \og(j) is of an even degree. 

Theorem 3.5. Let <p : 7L — > C be a finitely supported function such that there exists 9 q G 
/ = (—ir, tt], q G {1, . . . , Q} such that 

\4>(0 q )\ = 1 and |0(0)| <lonJ\{9 1 9 Q }. 

Assume further that for each q G {1,. . . ,Q}, there exist a q G K, j q G C with Re(^y q ) > 0, 
and an even integer v q such that 

0(9 q + 9) = 4>(9 q )e ia " e -^ e,Jq(1+om) at 9 = 0. (3.11) 

Then, for each N , there is a constant Cjv such that for any x G Z ; n = 1, 2, . . . , we have 

Q -. / I _ \\- N 
l^'MIS^E^ + t^P) . (3.12) 



Remark. An immediate application of Theorem 3.5 is that there exists a constant C such 
that, for all n, 

5> (B) (*)l<c- 



That is, we recover the positive part of the stability theorem, Theorem 2.1 Indeed, estimates 



such as (3.12) are more or less implicit in the proof of Theorem 2.1 given in [28, 29). 



Proof. We need to introduce smooth non-negative cut-off functions ip q such that each tjj q 
vanishes outside an interval 9 q + I q = [9 q — e q ,9 q + rj q ], ip q = 1 on [9 q — e q /2,9 q + rj q /2], 
9 q i [^ — 36^/2, 9 q + Zrj q /2], q' ^ q, and Y^i V'q = 1 (i- e -i the ^'s form a partition of unity). 
Using these cut-off functions, we write 



(n) 



(x) 



l - jf e— 9 H9) n d9 = E ^ / e-* xd m n ip q (9)d9 



2tt 

Q 

E 



— e- ix9 i4>(9 q ) n 



E 



E 

l 





2ir 




>4(9 g r 




2vr 


e -ix 


9 ^(^ 9 ) n 



-i(x—a q n) 



- l e- a " e 4>{9 q + 9)] n i/> q (6 q + 0)d9 



i(x ~ a ^ e P q (9) n 4> q (9 q + 9)d9 



2im l l v 



. (x — aqn) 



1/vqr 



e ^ P q>n (6)i> q , n (9)d9 



19 



where 



P q {6) = 4>(9 q r 1 e'°"'\ v „ l 
P q>n {9) = PgiO/n^r, 

i> q ,n{0) = ^ q {9 q + 9/n l ' V «). 



+ 



By hypothesis, 



P q {9) = e 



-761^(1+0(1)) 



at 0. 



The function P q can be viewed as an entire function of z = u + iv £ C which is periodic 
in u and has at most exponential growth in \v\ for large \v\. Hence there are constants 
7i)72 £ (0,oo) such that 



\P q (z)\ < exp(-7m^ +j 2 v Uq ) on {z = u + iv : u £ (-3e q /2,3r] q /2)} 



(3.13) 



and 



\P q (z)\ < exp(72^ 9 ) on C. 



(3.14) 



A priori, the constants 71,72 depend on q but since there are only finitely many q, one can 
assume that they are the same for all q. 
Next, for each q £ {1, . . . , Q}, consider 



1 



. (x-a q n) „ 

e p q , n {e)^ q ,n{o)de. 



Observe that 



. {X - OL q nY . 

n l/Uq ' 



1 



21m 1 1 v i J n 1 / V il 
1 



. (x — otqTl) 
— I 1 



P g AG)^A d )dO 



(x-a q n) 



e di[p q , n (9)^ q , n (e)]de. (3.15) 



The iterated integration by parts performed to obtain the second equality does not produce 
boundary terms thanks to the cutoff function V^n- The following proposition is analogous 
to Proposition 3.2 and the proof is the same. 

Proposition 3.6. There are constants C\,a\ £ (0, oo) such that for any q,n = 0, 1, 2 ... , 

and 9 £ (-n l / v ie q ,n l / v tri q ), 

\d k e P q , n (0)\ < Cl^klk-^expi-a^"). 



Using (3.15) and Proposition 3.6, we obtain 

k 



(x — a q n) 



n l/u q 



Ck 



(3.16) 



The desired result follows. The necessity to separate the contributions of the different 9 q via 
the use of cutoff functions prevents us to obtain more precise upper bound in this case. □ 
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4 Carne's transmutation formula 



This section develops a version of Carne's transmutation formula to obtain universal long 
range upper bound on M£(x,y) where M is a (finite range) normal contraction acting on 
£ 2 (X, 7r). Here, X is a countable space, tt a positive measure, M acts on £ 2 (X,tt) by 

Mf(x) = Y,M(x,y)f(y), 

y 

k is a fixed integer and 

M k = I-(I-M) k = Q AP. 

If we think of I — M as a " Laplacian" , then M£ is the discrete semigroup of operators 
associated with the k-th power of this Laplacian, namely, (7 — M) k . 
Assume that X is equipped with a metric d such that 

d(x, y) > 1 M(x, y) = 0. (4.1) 

This is a finite range condition in terms of the metric d. If the matrix M is finite range 
in the sense that for each x there are finitely many y such that M(x, y) ^ then we can 
define the metric dM by 

d M (x,y) = mf{k:M k (x,y)^0}. 



This metric obviously has property (4.1 ) (we do not have to assume that c?m is finite for all 
x,y). 

Our goal is to obtain bounds showing that MJ}(x, y) is small when d(x, y) is large when 
compared to n l ^ 2k \ In the classical case where k = 1 and M is the transition matrix of a 
reversible Markov chain, the Carne-Varopoulos estimate states that 

|M>,,)|<2f44V /2 exp^ ^ 



vr(y)y V 2n 

This bound is remarkable for its generality and explicitness. To see it in action for Markov 
chains on Z d , see [3]. For an extension to non-reversible chains, see [19]. For a probabilistic 
interpretation, see |22j . 

Our result reads as follows. 

Theorem 4.1. Let (X,d) be a countable metric space equipped with a positive measure tt. 
Fix k and sq G (0,2~ 1+1 / fc ). There exist constants C,c £ (0,oo) (depending on k and so) 
such that, for any normal contraction M on £ 2 (X,tt) satisfying (4.1) and whose spectrum 
is contained in [a, 1] with a E [—1, 1) and 1 — a < 2sq, 

Vx,yeX,n = l,..., \MZ(x,y)\<C (^^' 2 ^{-c ( ,U r ' n " 1L 



ir(y) J \ ^n 1 /^) 
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Remark. Proving that an operator is normal and computing its spectrum is not an easy 



task, in general. Most application of Theorem 4.1 are likely to involve cases when M is a 



hermitian contraction. In this case the hypothesis that the spectrum is contained in [a, 1] 
is a very natural hypothesis since the spectrum is real and contained in [—1,1]. 

Proof. Following [6], consider the Chebyshev polynomials 

Q m (z) = \((z + (z 2 - I) 1 ' 2 )™ + (z- (z 2 - l) 1 / 2 )™) , m e Z. 

Each Q m is in fact a polynomial of degree \m\. For a £ [—1, 1], let 

'2z-\-a 



Qa,m( Z ) = Qk 



I- a 



These are the Chebyshev polynomials for the interval [a, 1] C [—1,1]. See [6} Theorem 2']. 
As explained in [6l Sect. 2], assuming that the spectrum of M on £ 2 (X, tt) is contained in 
[a, 1], it holds that 



Qa,m(M) : £ 2 (X, tt) — > £ 2 (X, tt) is a contraction 



and 



M n = ^P (A£ = m)Q o , m (M) 

m 

where X% is the simple random walk on Z with transition probabilities 

P(A n = ±l|X„_i) = (1 - o)/4, P(A n = 0|X n _i) = (1 + o)/2. 
Consider the measures (3 S on Z where (5 is the Bernoulli measure /3(=tl) = 1/2 and 

P s = (1 - s)5 + s0, sG [0,1]. 



(4.2) 



Further, set 
By definition, 



Also, by (4.2) 



Similarly, 



A s / = / * (S - p s ), ^ k = 5 Q - (So - Psf k . 
F (XZ = m) = pi n Xm), s = (l-a)/2. 

Ml = (!-(!- MfT = £ ^(m)Qa,m(M). 
(I - MfMl = £ Ai^(m)Q a>m (M). 



(4.3) 



(4.4) 



Now, let (f>i, % = 1, 2, be two functions on A with support in = B(xi,ri) and assume that 
d(Bx,B 2 ) > r. Then 



((I-MYm^u^U = ^ { $(m)(Qa, m (M)<i)iA2).. 
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Further, since Q a ,m(M) is an £ 2 (X,tt) contraction and a polynomial in M of degree \m\, 

|(Qo,m(M)^i,^2) ff | < H01II2II02II2 

and 

(Qa,m(AO01>^2>w = if H < r. 

This last property follows from the hypothesis that M(x, y) = if d(x, y) > 1 which implies 
M l (x,y) = is d(x,y) > i>0. Putting these properties together, we obtain 

- mYM'Z^^A <\\<Pih\\foh £ lAV^MI, s = (l-a)/2. 

|m|>r 

For < s = (1 — a)/2 < so < 2~ 1+1 / fc , the Fourier transform 

&,fc(0) = l-8*(l-COS0)*) 



satisfies |i/! s &| < 1 on 7* and (3.3) with u = 2k and 7 = (s/2) k . One can conclude that 
there are constants C, c E (0, 00) depending only on s such that 

W-M)'M trh MA < ex P (-.(j^)*) ■ (4.5) 



The inequality stated in Theorem |4. 1| follows by taking </>i = 5 X , 4>2 = d~ X2 , x\ = x,X2 = y, 
ri, T2 = 0, r = y). □ 



Inequality ( |4.5[ ) contains more information than captured in Theorem 4.1 In particular, 
we get the following result. 

Theorem 4.2. Let (X,d) be a countable metric space equipped with a positive measure n. 
Fix k and sq € (0, 2~ 1+1 / fc ). There exist constants C,c £ (0, 00) (depending on k and so) 
such that, for any normal contraction M on £ 2 (X,tt) satisfying (4.1) and whose spectrum 
is contained in [a, 1] with a 6 [—1, 1) and 1 — a < 2sq, we have 



\ ( I-M)'K^v)\<C'[^j (l+n)-"exp „ a 



Further, if A x (r,R) = {y : r < d(x,y) < R} then 

C*7r(*) 1/2 

yi - ivi) ivi k yx,y)\ n^yj \ 

yeA x (r,R) 



E \V-M)'Mt( X ,y)Uy)< ^S»p(-<^t) 2 "') 



Although these results are far from optimal in the sense that they do not capture the 
decay of M'j}(x, x), they do contain some useful information as demonstrated in the following 
corollary which provide a highly non-trivial lower bound on M| ra (x, x). This follows closely 
one of the result obtained in [7] in the reversible Markov case and with k = 1. 
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Theorem 4.3. Let (X, d) be a countable metric space equipped with a positive me asur e 
it. Fix k and sq £ (0, 2~ 1+1 / fc ). Let M be a hermitian contraction satisfying (4.1), 



^2 y M(x,y) = 1 and with spectrum contained in (—2sq + 1,1]- Fix x £ X, set V(x,t) 
7r(B(x,t)) and assume that there exists a positive increasing function v such that 



V(x, 2t) < w(x)v(t), u(0) > (1 + 2C) 1/2 , 



and with 



1 1 — y t 2k - 1 1 log v(t) increasing to infinity. 
Given n, let r(n) be the smallest integer such that 

2k 



Then 



5 n V(2fc i) < cr 3k ~ 1 /log v(r). 

4 V(x, r(n)) 



In particular, 

1. IfV(x,t)<Air(x)(l + t) D then 

M% n {x,x) > avr(x)[(l + n)log(l + n)]- D / 2fc . 

2. IfV{x,t) < n(x)exp(At p ) with p e (0, 2k ^ 



'' 2k-U> 



M£ n (x,x) > ir(x)exp -A n 



g(2fc-l) 
,2fe(l-/3)+/3 



(4.6) 
(4.7) 

(4.8) 



Remark. The last inequality is informative only when f3 < l/(2k — 1) because M£(x,x) > 
ir(x)V(x, n) _1 , always. 

Proof. Set B(x, r) = {z : d(x, z) < r} and m^(x, y) = MJ}(x, y)/ir(y). Since M is hermitian, 
that is, M(x,y)/ir(y) = M(y,x)/n(x), and ~^2 y M(x,y) = 1, we have (the first step uses 
that M is hermitian) 

m 2 k n (x,x) = Y,\ m k(x,z)\ 2 Tr(z) > \ m k(x,z)\ 2 ^(z) 

z z£B(x,r) 

2 

rn%(x,z)n(z) 



> 



> 



ir(B(x, r)) 
1 

tt(B(x, r)) 



1- Yl l m ft(^)KO) 

z£B(x,r) c 
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Next, set A q = A(r2 q ,r2 q+1 ) and, using Theorem 4. 2l write 



53 K(*,*)K*) = ££ \ M k^ z )\ 

z&B(x,r) c 9=0 zeA, 

2k — 1 

1/0 „ (2 g r) 2fc-l x r(B(3; ,2';+ 1 1 -)) 1 

< CTr^^^e c „i/(2^-i) + 2 logl 1 



Using the function v given by (4.6)-(4.7) and any r such that 

5 n V(2fc) < cr sl^T / log w(r), 
the first term in the series above is bounded by 



2k 



e - c ^W+l 1 °g[ Z ^gr 11 ] < e -(9/2)log«(r)_ 
Further the ratio of two consecutive terms in the series is bounded by 

2k 

^l^jff^H^v^r) < e -21og,(r)_ 

Hence, 

C » e -(9/2)log«(r)+21ogu(r) 



e 21ogi;(r) _ ]_ 
z£B(x,r) c 



C C 1 

- u ( r )i/2( 7 ;( r )2 _ !) - u ( r )i/2( u ( )2 _ i) - 2" 



2n/ 1 7r(aj) 



It follows that 

Ml n {x,x) > A r> 

k y ' ~ 4ir(B{x,r(n)) 

where r(n) is the smallest integer such that 

5 n V(2fc) < cr sl^T / log w(r). 



□ 



5 Applications and history 

Iterations of signed kernels are useful in data smoothing, image analysis, density estima- 
tion, and summability. There are a variety of applications to partial differential equations 
via divided difference schemes for numerical computation and asymptotics of higher-order 
equations. These applications served as motivation behind most of the previously developed 
theory. In this section, we briefly give pointers to the literature and describe the basic ideas 
behind these applications. 
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5.1 Data smoothing, de Forest's problem, and density estimation 

If yi, ... , y n , ... is a sequence of numbers, it is standard practice to smooth this sequence 
by local averaging to eliminate noise. An early instance of this, due to C.S. Pierce, replaced 
Vi by y'i = {yi + ?/i+i)/2. This was iterated four times giving yf ] = 2" 4 J2*j=i (f)yi+j. 
Because the binomial distribution approximates a normal distribution, this is an early 
version of smoothing with a Gaussian kernel. This story is told by Stigler [26] . The next 
few paragraphs constitute a review of some more sophisticated examples. 



de Forest's problem 

It is natural to insist that a smoothing method applied to a constant sequence returns the 
same constant sequence. For convolution smoothers 



Vi 



(applied to two-sided sequences) this forces ^ f(i) = 1. In early work on smoothers, de 
Forest asked about smoothers that would preserve low order polynomials. To achieve this, 
one must allow / to change sign. For example, convolution by / supported on {0, ±1, ±2} 
with /(±2) = —1/9, /(±1) = 4/9 and /(0) = 1/3 preserves all quadratic polynomials, 
de Forest's problem refers to the study of the behavior of the iteration of such smoothers. 



See (5.4) below for this example. A charming history of de Forest with good references is 
in [25] and hard to find papers are reprinted in [26] . Schoenberg [23] studies de Forest's 
problem melding it with a discussion of variation diminishing and total positivity. This 
work is followed up by Greville \15\ [T6] . 



Other smoothing techniques 

In recent years, a host of techniques using iteration and negative weights have emerged. An 
early iterative technique called "3RSSH" smooths a series by replacing yi by the median of 
yi-l, yi, yi+i [30]. This is iterated until it is stable. Tukey introduced a smoothing method 
called "twicing" (see, e.g., [U Sect. 11] and \12\ I27j). This is a linear method in which a 
series yi is smoothed with a kernel K and then the residuals yi — y[ are smoothed by K 
and added in. This amounts to smoothing with 2K — K * K. A further iteration leads to 
3K — 3K *K + K*K*K. More generally, I - (I - K)* k gives the fc-fold iterate. This fits 
exactly into some of questions discussed in Proposition |2.4| Indeed, Tukey's smoothing is 



based on the Bernoulli measure (3 (with no holding). Proposition 2.4 shows that the n-th 



convolution power of = 2/3 — (5^ blows up with n. The Bernoulli measure f5\ii (which 
has holding equals to 1/2) provides a measure such that the convolution powers <pi of 
4>k = <5o — (£o — ^\j2) k behaves well for all k. 

Density estimation 

Given a random sample y±, . . . , y n from an unknown probability density g, the kernel density 
estimator of g is 
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A huge statistical literature studies the properties of such estimators, their behavior when 
n is large, choice of h and choice of the kernel K. A standard choice is the Gaussian kernel 
K(x) = ~^= e ~ x ^ '■ A recurring theme known variously as "higher-order kernels" or "super 
kernels" allows kernels that change sign. A kernel is of order k if 



/ 



x m K(x)dx 



1 if m = 1 

if m = 1, . . . , A; — 1 

c if m = k (c^O). 



Equivalently (assuming K has high moments), the Fourier transform K satisfies K(£) = 
l + c£ fc + o(|£| fc ). The paper [T7] develops these ideas in modern language and focuses on 

K{x) = (27T)- 1 J cos(te)e" l * |fe dt, K{0 = e'^" , k > 2. 

For even k, these are exactly the de Forest-Schoenberg-Greville limiting kernels. Many 
variations of these ideas have been explored and [U [13] are good sources for this material. 
The associated density estimators have better mean-square error but, because higher order 
kernels must have negative values, may give negative estimators for g. 

5.2 Partial differential equations 
Higher-order evolution equations 

Natural scientific problems give rise to higher-order equations such as the initial value 
problem of solving for u(x,t) satisfying 

du 

— = Hu, given u(x, 0) = v(x) for x £ M. n , 

with H, say, a uniformly elliptic operator of order 2m densely defined on L 2 (M n ). Under 
conditions on H, u(t, x) can be expressed as 



i(t,x) = (e m v)(x) = [ K(t,x,y)v(y) 



dy. 



If H has constant coefficients, K(t,x,y) = P t (x — y) = j^y; J Rn e l ( x y ^ d£ with the 
polynomial the symbol of H. In the special case that p(£) is homogeneous of order 
2m (so p(a£) = a 2m p(£)), P t (x) = t- n / 2m P{xt- l l 2m ) with P = P x . Thus estimates on 
P{ x ) = (2^)" Ir" e lx £~ p ^ dl; determine the long-term behavior of K and so u(t,x). This 
program is explained and carried out in a series of papers by Brian Davies and coauthors 

E3 ei Eg maun. 

Consider the case where n = 1 and Hu = ^s%u. Then p(£) = £ 2k and P{x) = 
5¥ Jk e lx ^~^ 2k d£. This is exactly Pi2k{x) of (1.3). It is well known that 



H ( ,_ 1 y(-^) 2 r(2n + l) (Q , T (^) r H M-1 
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Further there exists a,b,c E (0, oo) such that for all m = 1, 2, 3, . . . and real u, 



d 

du 



H 2 k{u) 



< c 



rn^ m e _ aM 2fc/(2fe-u 



See |144 Chap. 4]. This is also true (with 0° = 1) with m = 0. Pictures of H2 and H4 are 
in |Figure 1 of Section 1 



Divided difference schemes 

Iterated convolutions of complex kernels arise in finite difference approximations to the 
solutions of partial differential equations. Consider first the simple equation 



du du 
dt dx J 



0<t<T, < x < 1, u(x, 0) = v(x) given. 



(5.1) 



In (5.1) a is a given constant and u{x,t) is to be found given v. A finite difference scheme 
for this equation discretizes time into multiples of At and space into multiples of Ax. If 
Xi = iAx, tj = jAt, one simple scheme is 



u(xj,t j+ i) - u(xj,tj) _ a(u(xi +1 ,tj) - u(xj,tj)) 
At Ax 

or, with A = At/ Ax, 

u(xi,tj + i) = (1 - a\)u(xi, tj) + a\u(x i+ i,tj) = <jP> * v(xi) 



(5.2) 



(5.3) 



with 0(0) = (1 — aX), tf>(— 1) = aX. Thus, "running" the differential equation corresponds 
to repeated convolution of the (perhaps complex) probability measure (p. Higher-order con- 
stant coefficient equations and more sophisticated difference schemes lead to more complex 
functions (j). A readable introduction to finite difference equations is in |23| . 

One may approximate even-order derivatives by symmetric differences; writing h = Ax, 



fix - h) f(x + h) 

o f{x) o 



f (2m) (x) 



1 

¥ 



{f{x - 2h) - 4f(x -h) + 6 f{x) - 4/(x + h) + f(x + h)} 



i=o v J / 



(j -m)h). 



Similarly, when ^(1) = — 1) = \, /u(x) = otherwise, we have 



x 


-1 





1 


6 Q - n 


1 
2 


1 


1 

2 



and 



Thus, 

f(2) = 



X 


-2 


-1 





1 


2 


iSo - v) {2) 


1 

4 


-1 


3 
2 


-1 


1 

4 



j^ii -n)*f, / (4) = %ii-») {2) */.-•-, f [2m) 



(-2Y 



-iSo-v) {m) *f. 
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Now consider the equation 

du d 2m 



dt dx 2m 

Discretizing, 



u with u(x, 0) = v(x). (5-4) 



or 



u(xj,t j+1 ) -u(xi,tj) = (~2) m _ (m) ^ 
At /i m 1 W 

ufc, ij'+i) = (<5 + A(5 - /x) (m) ) * «(xi, ii), A = (-2) m At/h r ' 



If ft, = 2(At) 1 ^ m , this gives again 5o — (Sq — /i)( m ) = <fi, thus again, running the differential 
equation ( |5.4| ) corresponds to repeated convolution by the signed probability 4>. 

A more sophisticated discussion of these ideas together with assorted results can be 
found in Thomee's survey 
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